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NAVAL ARCHITECTURE 


The catamarans George and Energa Solar 


Dariusz Duda, 
Wojciech Litwin 
Gdansk University of Technology 


ABSTRACT 


Since a few years students of the Faculty of Ocean Engineering and Ship Technology, Gdansk University 

of Technology , have designed and built untypical floating units. Until last year their efforts were focused 

mainly on leg-driven boats. The boats of interesting design have taken part in yearly competitions : the 

International Waterbike Regatta. Their advanced design made it possible to compete with the best boats 

from Croatia, Holland, Germany, Turkey and Italy. Recently the students have designed and built a solar- 

energy-driven boat. It is the catamaran Energa Solar which took part in a prestigious regatta : the Frisian 
Nuon Solar Challenge carried out in Holland in summer 20006. 


Keywords: leg-driven boats, solar- energy - driven boats, renewable energy 


INTRODUCTION 


The Scientific Circle of Students Korab has acted for many 
years at the Faculty of Ocean Engineering and Ship Technology, 
Gdansk University of Technology, however presently an extra- 
ordinary activity was shown by the students associated in it. 

The revival happened already in 2004 when a group of stu- 
dents attempted to design and build a leg-driven boat intended 
for taking part in the International Water-bike Regatta whose 
6" edition had to be performed in Bremen in May 2005. 

The RW-4 catamaran designed and built within a short 
time, made successively its debut : it was ranked sixth among 
28 boats built by students from Croatia, Holland, Germany, 
Turkey and Italy. The success has motivated the team to an 
even greater effort. Basing on experience gained during the 
regatta it was decided to build an entirely new boat which could 
be able to compete for a place on podium. In the end of 2005 
the Rector of Gdansk University of Technology obtained an 
official inquiry from the side of the organizers of the regattas 
of solar bikes, Frisian Nuon Solar Challenge, if somebody of 
the University would be interested in starting. The students of 
the Scientific Circle promptly declared their will to take part. 
However, as revealed clear in a short time, the task imposed 
on themselves, namely the building of two boats practically in 
parallel appeared difficult to be obtained, as the International 
Waterbike Regatta was scheduled for the end of May and the 
Frisian Nuon Solar Challenge — for the end of June, i.e. during 
examination session. 


DESIGNING AND BUILDING THE 
CATAMARAN GEORGE 


Debuting in 2005, the catamaran RW-4 was reliable and 
fast, and of excellent manoeuvrability due to its azimuthing 
propeller. Its main drawback was large weight resulting from 
the application of simple polyester composites (for hulls) and 
usual structural steel (for frame). 


In designing the new boat, was taken into account the spe- 
cificity of the regatta within the frame of which seven highly 
different competitions such as : sprint, slalom, acceleration 
test, forward-and-back sprint, long distance run, bollard pull 
test, are carried out. 

Therefore mass reduction of the boat and achieving its 
possibly high speed even in expense of lowering its manoeu- 
vrability were assumed design priorities. 

The analysis of the competed structures revealed that each 
of them had important disadvantages. Therefore the team pre- 
sented three different novel design concepts. It was decided 
to build a catamaran driven by a pulling azimuthing propeller 
placed fore. To limit number of mechanical gears and to unload 
the boat structure the crew members were assumed to occupy 
places along the hulls (Fig. 1). 


Fig.1. The catamaran George . 


The team planned to single-handedly design and build 
the hulls but cost and time consumption of such undertaking 
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forced it to drop out the plan. The model tests realized by Ship 
Design and Research Centre (the CTO) in Gdansk (Fig. 2), 
demonstrated that the previously applied hulls of typical canoe 
form are of a low resistance at an assumed speed of motion. 
Hence it was decided to use a more modern canoe hull. Owing 
to the Gemini company which made available hull moulds the 
students unaidedly produced two laminated hulls. 


Fig. 2. Model tests carried out by the CTO. 


The designed structural frame was computationally con- 
trolled by using the Finite Element Method. Most troubles 
were associated with the highly loaded pipe connecting both 
hulls in the bow part of the catamaran. It transferred wave- 
-generated loads and torque and bending due to azimuthing 
propeller thrust. 

The connecting frame was built of a high-grade aluminium. 
All frame elements were carefully prepared by the students 
and welded together by the company Aluminium Ltd at Wista 
shipyard. 

The propulsion system was so designed as to reduce, as 
much as possible, torque in the azimuthing propeller vertical 
shaft. It is important since the torque interacts with the pro- 
peller rotation mechanism that may be troublesome for swain. 
For this reason the rotational speed of the first two stages of 
the reduction gear was increased up to the assumed rotational 
speed of the propeller. As demonstrated later during practical 
tests, the torque acting in the rotation mechanism was negligible 
and did not make any difficulty for swain during manoeuvres. 
The body of the intersecting-axis gear which operated under 
water, was modernized for two reasons. Its flange connection 
with the propeller column, provided for by its producer, was 
feasible but of large gabarites, that could detrimentally influ- 
ence resistance to motion of the immersed part of propulsion 
system. Hence the reduction gear was disassembled and the 
upper connecting flange removed out of the body by machi- 
ning. Then the propeller column was glued with the gear body. 
Additionally it was assessed that the propeller was to be moved 
away from the propeller column, hence the existing propeller 
shaft was replaced with a longer one and one sliding radial 
bearing was added. 

The propeller, a crucial element of the boat, was produced 
by the CTO whose help appeared invaluable. The catamaran’s 
design process took a few months, but its building and assem- 
bling — only four weeks. 

First trials performed on Motfawa river showed that the 
catamaran’s centre of gravity should be shifted aft as the boat 
was slightly trimmed by the head. 

After some modernization work resulting from the trials 
the catamaran was finally completed one day before leaving 
for the regatta. 

After traveling via Bohemia, Slovakia, Hungary, Serbia 
and Bulgaria the team finally reached Turkey. There after hard 
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Fig. 3. The simplified assembling drawing of the azimuthing propeller with 
vertical shaft; Variant A — for the catamaran George, Variant B — for the 
Energa Solar; 1 — seating of additional bearing supporting the shaft, 

2 — the longer shaft made of stainless steel, 3 — additional sliding bearing 
supporting the shaft, 4 — modernized intersecting-axis reduction gear, 

5 — vertical propulsion shaft, 6 — propeller column fitted with bearings 
and rotation mechanism flange, 7- rotation mechanism body, 

8, 9 — propulsion engine. 


fighting it managed to win 1* place in general classification. 
The challenge cup founded 28 years ago by the German team, 
first time reached Poland. 

The team will face hard work in the next year as the Inter- 
national Waterbike Regatta of 2007 has to be held in Gdansk. 


DESIGNING AND BUILDING THE 
CATAMARAN ENERGA SOLAR 


To be carefully acquainted with the rules was the first 
step in preparation of the team, made just after submitting its 
participation in the regatta. The regulations of over twenty 
pages were very extensive and detailed. It imposed values of 
maximum dimensions of the boat, its mass and number of crew 
members, as well as it even determined the minimum weight 
(mass) of each crew member amounting to 70 kg. 


The regatta organizers established 
three separate classes of boats : 


* The Class A — for one-person boats having solar cells of 
less than 100 kg mass and 6 m length at most 

* The Class B — for two-person boats having solar cells of 
less than 150 kg mass and 8 m length at most 

* The Open Class —for boats of 8 m length at most, without 
any other limitations. 


As Dutch partners suggested that the Class A would have 
the greatest number of competitors it was decided to build an 


one-person boat complying with relevant requirements for 
that class. 

The team faced difficulties in gaining financial support for 
building the boat. As in the case of the catamaran George it 
was the Rector of Gdansk University of Technology and the 
Head of the Faculty of Ocean Engineering and Ship Technology 
who financially supported the building of the boat. However it 
was crucial to find a strategic sponsor. Ultimately, the company 
Energa S.A. accepted position of the main sponsor. 

To take part in the regatta it was necessary to fulfill a few 
conditions. The regatta’s organizer required a.o. to submit the 
boat’s technical plans and schematic diagram of electric ne- 
twork. The team was concerned about if crew safety is the only 
thing of this request or perhaps also the gaining of experience 
on the basis of the submitted designs. 

It was important that the company Sharp declared to equip 
each team with the same number of identical solar cells. They 
were only lent and had to be returned after the regatta. However 
this way all the teams had the identical energy sources, that 
- in the conditions of limited funds - was crucial and helpful in 
getting even to each other. The energy source for the A-class 
boats consisted of six panels whose total power output could 
be up to 600 W depending on solar operation. 

The regatta’s organizers allowed to store the produced 
energy in electric batteries of the capacity limited to about 40 
Ah at 24V voltage. No charging them from land was allowed 
and only solar radiation had to be used. 

The organizers attached high importance to safety issues 
therefore each boat had to be fitted with a fire extinguisher, water 
draining pump, emergency electric ,,death man switch” etc. 

The planned route of the regatta run along canals of Holland. 
It was divided into six stages and its total length reached 230 
km. In designing the boat it had to be remembered that the boat 
must be durable and reliable enough to cope with several day 
trip in variable ambient conditions. 

As in the case of the catamaran George, the boat’s concept 
was discussed for a relatively long time. It was sure that an 
azimuthing propeller will be applied. The team has already 
gained some experience concerning the boats of similar mass 
and propulsion power. Hence, knowing advantages of the 
catamaran concept and necessity of fitting six cell panels of a 
mass close to that of the boat itself, the team was convinced 
that the twin-hull unit should be chosen. The problem was to 
get appropriate hulls as the firm which already accepted the 
building order withdrew from it in the last moment. The hull 
of the canoe used up to that moment appeared to short and its 
displacement insufficient. In this situation 6 m hulls of one of 
typical sailing catamarans was chosen. The order for two such 
GRP hulls was placed with one of their producers. To limit 
their weight thickness of laminate was kept at minimum and 
original heavy decks were removed. As prompt delivery of the 
hulls was not possible it was also not possible to finally design 
and build the frame, without which to determine location of 
centre of buoyancy necessary for arranging solar cell panels 
and position of swain, was impossible. 

Two months before the regatta the delay connected with 
building the frame and hulls, increased since it was necessary 
to focus the team’s effort on the water bike because the star- 
ting date of the International Waterbike Regatta was nearer 
and nearer. 

As late as in the end of May the work upon the solar - energy 
- driven catamaran started again when only four weeks rema- 
ined to the date of departure for taking part in the regatta. The 
approaching examination session made the situation even worse 
as the students had less and less time to build the boat. Never- 
theless the frame and hulls were completed on schedule. 


An azimuthing propeller was used to drive the catamaran. 
From the very beginning of designing process of the boat two 
different design solutions were considered (Fig. 1 and 2). It was 
possible either to place the electric motor within an underwa- 
ter pod and apply direct propulsion or to drive the propeller 
through an intersecting-axis reduction gear and vertical shaft 
and to place the motor over the propeller’s column. Tests of 
the electric motor on a special test stand showed that to apply 
a rotational speed reduction gear would be more favorable. 
Unfortunately the directly driven azimuthing propeller was al- 
ready ready to use and its modernization was rather impossible. 
Time was running and funds lacking. Therefore elements of 
the propeller with vertical shaft was promptly ordered and an 
intersecting -axis gear able to reduce rotational speed in half, 
was purchased. 


l 2 3 4 6 


Fig. 4. The simplified assembling drawing of the azimuthing propeller with 
the electric motor placed in the pod; 1 — fairwater, 2 — motor 5 seating, 

3 — driving motor, 4 — stainless steel shaft with pair of bearings, 5 — second 

part of the propeller s pod with the seating of rolling bearings, 6 - propeller . 


To reduce the propeller’s resistance to motion the reduction 
gear was subjected to the modernizations already proved in the 
propulsion system of the catamaran George (Fig. 1), i.e. the 
flange connecting the pod with propeller’s column was removed 
and the propeller itself shifted away from the reduction gear. 


Fig. 5. The catamran Energa Solar . 


To drive the propeller a non-typical three-phase electric 
motor fitted with permanent magnets containing neodymium, 
was used. Owing to its relatively large diameter and a large 
density of magnetic field generated by the permanent magnets 
it produced a relatively large torque. In spite of its small mass 
of 680 g its output power reached almost 2k W. Results of tests 
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performed on a special test stand showed that the motor oper- 
ated correctly within the assumed range of its rotational speed. 
It was assumed that during the regatta the motor has to oper- 
ate on the power level of 500 — 700 W. Unfortunately, during 
the tests some troubles associated with the motor’s controller 
overheating, appeared. The controller fed from DC source 
(accumulator battery) produced three-phase alternate current 
feeding the motor. Its efficiency was relatively high reaching 
even about 90%. In effect large amount of heat was emitted. 
Due to compact design of the controller, heat exchange was dif- 
ficult and resulting increase of temperature triggered the thermal 
cut-out to operate until the temperature drops enough. Because 
of lack of time it was not possible to order another controller. 
Hence the previously purchased controller was subjected to 
modernization work. The controller’s circuit plate packets were 
divided into three separate segments and glued onto a radiator 
fitted with a fan (Fig. 6). The solution appeared effective and 
turned out to be correct during the regatta. 


Fig. 6. The DC/AC controller; on the left, inside casing - its modernized 
version, on the right — its original version before modernization . 


One of the prerequisites for taking part in the regatta was 
to submit to the organizers the electric network documentation 
for careful review in advance the regatta. Necessity of such 
control has been demonstrated by fires broken out on two boats 
as aresult of overloading the electric network. The electric 
network of the catamaran was based on two 12V gel batteries 
connected in series. Each of them was placed in another hull. 
Left from the swain post a box containing all electric and 
electronic devices, was located (Fig. 3), except of the motor’s 
controller which was placed near the driving motor. Initially 
it was planned that in the case of propulsion system’s failure 
the fast switching from 24 V to 12 V voltage in the electric 
network, would be possible. Such option could make it po- 
ssible to install a typical overboard motor used for propelling 
small boats. However tests carried out in the model basin of 
the CTO demonstrated that the catamaran’s speed fitted with 
a 300W outboard motor would be very low, hence the concept 
has been abandoned. 

To drive the boat a special propeller designed and manu- 
factured by the CTO was applied. 

In the course of the catamaran’s manufacturing efforts 
were undertaken to reduce its weight as much as possible. In 
spite of that its mass was exceeded by almost 15 kg. It was 
expected that the regatta organizers would penalize the team 
with penalty minutes, but they came to the conclusion that the 
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Fig. 7. Schematic diagram of the electric network of the catamaran Energa 

Solar; 1 — module of AC/DC controller; 2 — controller charging accumula- 

tors, 3 — cut-out in solar cells circuit, 4 — cut-out in circuit of accumulators, 

5 — cut-out in feeding circuit, 6 — power measurement system of solar cells, 
7 — power measurement system of motor, 8 — safety cut-out switch . 


overweight would make harm to the team itself and no penalty 
was finally given. 

A week before the regatta the boat was two times practically 
tested. Results of the tests appeared positive : the boat was fast 
and responsive. It was given the name of Energa Solar. 

The preparations were continued till the day of departure. 
Despite some troubles the team succeeded in being in time 
for the regatta. The boat, after thorough technical control, was 
approved for taking part in the regatta. 

Six days of the race were full of dramatic events. What’s 
most important the team succeeded in finishing the race in 
spite of some technical problems. Unfortunately the organizers 
punished the team for introducing the changes to the electric 
system. The penalty wrote off chance of winning a place on 
podium and finally the team was ranked seventh among 14 
boats of A class. 
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Application of fuzzy inference to assessment 
of degree of hazard to ship power plant operator 


Tomasz Kowalewski 
Antoni Podsiadto 
Wiestaw Taretko 

Gdynia Maritime University 


ABSTRACT 


This paper presents application of fuzzy logic to assessment of degree of hazard to ship power plant ope- 

rator. For the assessment a system of computer-aided identification of hazardous zone within ship power 

plant, was used. The system 5 variables representing the subject-matter knowledge in safety design area 

were transformed into fuzzy sets by means of appropriate linguistic variables and membership functions. 

The assessing of safety level of operator with the use of fuzzy inference was performed by means of an 
expert system programmed in the PROLOG LPA language. 


Keywords: safety, operator, ship power plant, assessment, fuzzy inference 


INTRODUCTION 


Complex technical objects (e.g. large transport units, 
electric power plants and power systems, chemical systems) 
can create significant hazards to safety of their operators. Sea- 
-going ship is a complex technical object specific from safety 
point of view. According to many domestic and worldwide 
elaborations and reports the ship power plant can be deemed 
the most hazardous place on ship. 

In ship power plant many dangerous and noxious factors 
which constitute various hazards to its operators appear simul- 
taneously. For instance the following can be listed: displacing 
machines and transported objects, moveable elements, falling 
elements, pressurized fluids, slippery and uneven surfaces, 
limited spaces, hot and cold surfaces, caustic and toxic substan- 
ces. Therefore it would be advisable to indicate, just in early 
design stages, such places in a ship power plant, which create 
potential hazards to operators. 

The question appears whether it would be possible to limit 
consideration of ship power plant only to potentially most 
hazardous zones during designing its safety. In the authors’ 
opinion it is possible provided use will be made of informa- 
tion on structure and function of a power plant, contained in 
preliminary design documentation, as well as its consideration 
will be limited only to the zones where an operator (operators) 
carrying out certain service operations concerning a given 
technical object, appears. The idea was realized in the computer 
aided advising system for identification hazardous zones in ship 
power plant [1]. The system proposes certain safety design 
strategies for the zones distinguished by it. 


The system’s concept consists in cooperation 
between designer and computer where : 


@ the designer provides appropriate information to the sys- 
tem on the basis of analysis of preliminary design of ship 
power plant, as well as his own knowledge, intuition and 

experience 


@ the computer processes the put-in data, calculates indices 
and orders distinguished hazardous zones of the power plant 
taking into account possibility of creating potential hazard 

to its operator. 


The description of the computer-aided identification system 
of hazardous zones within ship power plant was presented in [2], 
its elements (variables and functions) as well as a way of model- 
ling the subject-matter knowledge for computer purposes - in [3], 
and a way of determining its decision variables - in [4] , [5]. 

In the identification system in question the additive method 
of estimation of degree of hazard to ship power plant operator, 
was chosen. A drawback of the method is the linear ranking of 
potential influence of particular states (symptoms) of distin- 
guished variables, which could not be capable of modelling 
real influence of endangering and noxious factors on creating 
a hazard to operator. The drawback can be eliminated by using 
the fuzzy assessment and inference method. It correctly models 
the uncertainty associated with the assessment of influence of 
particular variables on the creating of a hazard to operator. 

This work is aimed at consideration of possible application 
of fuzzy logic to assessing hazardous zones to ship power plant 

operators. 


THE COMPUTER - AIDED 
IDENTIFICATION SYSTEM 
OF HAZARDOUS ZONES 
WITHIN SHIP POWER PLANT 


Operation of the computer - aided hazardous 
zone identification system 


The operation of the system in question is realized in such 
a way that its user — making use of a knowledge base and 
taking into account various limitations — performs assessment 
of the system’s state, and on this basis he determines a degree 
of hazard to operator by using a set of measures. The decision 
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—making procedure consists in choosing the values of input 
variables, which are capable of causing a definite effect to 
operator, i.e. hazard to his life or health. Basing on the values 
of input variables — determined by the system’s user — the 
system calculates value of hazard degree assessment index for 
considered element of hazardous zone and for n-th identification 
task, in accordance with the following relationship : 


J 
Mms Datus Kn, mWnCn + WjCjm (1) 
j=1 
where : 

a — index of influence of endangering and noxious factors 
on the creating of potential hazard to operator, which 
occur during carrying out service operations realized 
in compliance with the considered n-th procedure 
concerning m-th structural unit 

— index of influence of service factors on the creating 
of potential hazard to operator, which occur during 
carrying out service operations realized in compliance 
with the considered n-th procedure concerning m-th 
structural unit 

— index of influence of functional factors on the creating 
of potential hazard to operator, which occur during 
carrying out service operations realized in compliance 
with the considered n-th procedure concerning m-th 
structural unit 

— number of service operations appearing in n-th pro- 
cedure for m-th structural unit 

— weighing factor for n-th procedure 

value of service variable for n-th procedure 

— weighing factor for j- th functional variable 

value of functional variable for m-th structural unit 

number of functional variables. 


(n) 
LoF,m 


le ,m 


If the task of identification of hazardous zone element is 
considered for a set of service procedures it will be necessary 
to calculate the hazard index in compliance with the follow- 


ing realtion : N 


= (n) 
lonm= >. OE m (2) 
n=l 
where : 
I — total index of influence of service factors on the 


OF,m n + r: 
creating of potential hazard to operator, which occur 


during carrying out all service operations distinguished 
by considered procedures concerning m-th structural 
unit 

N — number of considered procedures. 


The assessment index of hazardous zone element — con- 
sidered according to a set-in criterion (number and kind of 
service procedures) — is determined by using the following 
relationship : 


ILogm Ikem 
lym= y +- (3) 
b2 Torm > Term 
m=1 m=l 
where : 


M — number of structural units distinguished 
by the considered procedures. 


Decision variables of the hazardous zone 
identification system 


As assumed according to [1], ship power plant operator 
will find himself in a potentially hazardous zone only when 
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he performs definite service operations. Degree of hazard is 
a function of factors resulting from operation of machines and 
devices, access to working place, kind of performed operation 
as well as position of operator. The first two are associated 
with functional structure of power plant and constructional 
form of considered unit and its environment, whereas the two 
remaining — with kind of service task carried out in given 
external conditions. For this reason the sets of the factors were 
conventionally split into the group of functional factors and that 
of service factors. For purposes of the identification system in 
question both the groups of factors were transformed into the 
set of input variables for the system. 

The performed classification of the factors hazardous and 
noxious to operator as well as the awareness of main and au- 
xiliary working processes carried out in ship power plant, as 
well as of service states of ship and its power plant, made it 
possible to distinguish the set of 10 input variables to the ha- 
zardous zone identification system, where 6 of them belong to 
the group of functional factors, and 4 to that of service factors. 
The detail description of the procedure for distinguishing the 
input variables is given in [4] and [5]. 

In this paper only the group of functional factors is consi- 
dered. The names and states (symptoms) of functional input 
variables used for assessing particular structural units are 
presented in Fig.1. 

In the system in question to each of the distinguished 
variables is attributed the weighing factor w, of constant value 
and five states (symptoms) of the values : 1, 0.75, 0.5, 0.25, 
respectively, or 0 in the case if a given symptom of considered 
variable does not involve a hazard to operator. 

The task ofthe system’s user is to attribute — on the basis of 
his knowledge, experience and intuition — appropriate values to 
the variables during solving a chosen identification problem of 
zone hazardous to operator. In the case when a variable takes 
zero value its symptoms are not taken into account on the sys- 
tem’s dialogue screens (Fig. 1) and then the system’s user does 
not mark any of those symptoms whereas the system automati- 
cally attributes zero value to the variable in question. 


FUZZY INFERENCE ASSESSMENT 
OF DEGREE OF HAZARD 
TO TECHNICAL OBJECT’S OPERATOR 


In considering hazard to operator all influencing factors, 
namely functional and service ones should be taken into acco- 
unt. However for purposes of analysis of possible application 
of fuzzy inference to assessing degree of hazard to operator 
the authors’ attention was paid only to functional factors, 1.e. 
those associated with realization of operational processes in 
various service states. To this end, the functional variables 
were transformed into linguistic ones and fuzzy estimates were 
attributed to their particular symptoms. Next, to particular 
linguistic variables appropriate membership functions were 
attributed. Trapezoid is one of the most often appearing forms 
of membership functions [7, 8, 9]. Hence the membership func- 
tions in question were represented by trapezoids. The assumed 
form and distribution of membership functions constitutes 
a preliminary proposal only. In the next phase of the research 
the properties of the functions will be corrected depending on 
the knowledge gained from the experienced experts — ship 
power plant operators. The particular membership functions 
for the decision variables (V1+V6) as well as the resulting 
variable which characterizes level of hazard to operator are 
presented in Fig.2. 

To assess degree of hazard to operator were used the pro- 
gramming language Prolog and software Flint which makes 


Assessment of the selected ship machinery component: “Transporting pump” 
identified in the following ship machinery items: “Fuel oil piping system — H.F.O. transfer piping system 
and involved in carrying out the selected maintenance procedure: Replacement of transfer pump 


” 


Heat energy load hazard — regular main engines, generator sets and auxiliary boilers 
operation of the machinery main engines, generator sets or auxiliary boilers 
component is necessary to assure main engines 
proper work of: generator sets or auxiliary boilers 


more than 3.0 MPa 
between 1.5 and 3.0 MPa 
between 1.0 and 1.5 MPa 
less than 1.0 MPa 


Fluid energy load hazard — the 
fluid pressure is: 


more than 200°C 


Thermal hazard — the fluid between 100 and 200°C 
temperature is: between 50 and 100°C 
less than 0°C 


exhaust gases 


o 

Chemical hazard — possible o fuel 

contact of operator with: o oil, cooling media, sanitary sewage 
o oiled, cooler or sea water, steam 


o participating in transformation of rotating motion into 
reciprocating motion 
Dynamic force hazard — possible o executing reciprocating motion 
contact of operator with elements: | o executing rotating reverse motion 
o executing rotating unidirectional motion 


converting operating/warning information and converting 
electrical energy into mechanical and thermal energy 
Electrical hazard — the electrical converting operating/warning information and converting 


ship machinery components are electrical energy into mechanical or thermal energy 
connected with: converting electrical energy into mechanical or thermal 

energy 

converting operating, warming or alarming information 


Fig. 1. Schematic diagram of assessing hazards to operator from the side of the structural unit : 
“transporting pump” during realization of the procedure: “replacement of transporting pump” . 


high 


Membership 
Membership 


Membership 
Membership 


cold dangerous 


Membership 


100 
temperature [°C] 


150 30 60 90 120 150 180 


electric current voltage (AC) [V] 


high 


hazard level 


Fig. 2. Membership functions for functional variables and hazard level . 
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fuzzy inference in that language possible. The base of rules in 
the software Flint, for the distinguished linguistic variables and 
their values, has the form of the fuzzy matrix as follows : 


fuzzy_matrix(t) :- 
V1 * V2 * V3 * V4 * V5 * V6 -> hazard level 


small*low*normal*small*small*safe -> low 


The notation can be interpreted as follows : „If V1 
is small and V2 is low and ... and V6 is safe, 
then the hazard level is low”. 


However in this case to fully determine all hazardous situ- 
ation it was necessary to write as many as 5120 rules. Owing to 
the great number of variables the base of rules was simplified in 


such a way that hazard level is considered separately for each + V1 (range 0+10) — 3 
of the variables, that significantly decreases size of the base + V2 (range 0+3.5) 31.25 


of rules. For instance the rules concerning the variable Z2 will 
have the following form : + V4 (range 0+10) > 6 


10 


fuzzy_matrix(pressure_reg) :- 
pressure -> hazard level 
low -> low During fuzzy inference process, due to aggregation of 
active rules, is obtained the resulting fuzzy set (Fig.4) from 
which — after its “defuzzification” — the final value of hazard 
v.high -> high. level is achieved. 


medium -> medium 
high -> significant 


+ V5 (range 0+10) > 2 
+ V6 (range 0+380) — 220. 


Assessment of the selected ship machinery component: “Transporting pump” 
identified in the following ship machinery items: “Fuel oil piping system — H.F.O. transfer piping system” 


and involved in carrying out the selected maintenance procedure: Replacement of transfer pump 


Heat energy load hazard (V1) — this factor concerns ship machinery components to be involved in 
converting chemical energy (contained in fuel) into mechanical, electrical or heat of energy, 
for instance: a fuel pump, an air cooler, pipes, etc. 


Estimation of hazard degree connected with participation in working process 
execution (in range 0 + 10): 3 


Fluid energy load hazard (V2) — this factor concerns ship machinery components to be involved in 
transporting high pressure fluids, for instance: a ballast water system, a bilge pump, etc. 


Estimation of fluid pressure (in range 0 + 3.5 MPa): 1.25 


Thermal hazard (V3) — this factor concerns ship machinery components to be involved in 
transporing heat fluids with temperature dangerous for operators, for instance: an air cooler, 
a exhaust pipe, fuel delivery pipes, etc. 


Estimation of fluid temperature (in range -50 + +250°C): 150 


Chemical hazard (V4) — this factor concerns ship machinery components to be involved in 
converting or transporting chemical substances dangerous for operators, for instance: pipes and 
pumps for transporting fuel, oils, etc. 


Estimation of chemical hazard degree (in range 0 + 10): 6 


Dynamical force hazard (W5) — this factor concerns ship machinery components whose elements 
execute reciprocating or rotating motion, for instance: engines, pumps, compressors, etc. 


Estimation of hazard degree resulted from dynamical force influence 
(in range 0 + 10): 


Electrical hazard (V6) — this factor concerns ship machinery components whose elements are 
charged by electrical energy, for example: switchboards, electric motors, etc. 


Estimation of electric current voltage - AC (in range 0 + 380 V): 220 


Fig. 3. Functional variables and estimates involved by the system Ss user . 
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Full estimate of hazard to operator by using fuzzy sets, 
performed by the system’s user for a given structural unit, will 
be based on the estimates of particular functional variables. 
The partial index of hazard level assessment for the considered 
hazardous zone element, I, ,,, will result from the fuzzy infe- 
rence. The value will be obtained by means of the MAX-MIN 
inference method with implemented ,,defuzzification” by using 

the gravity centre method [7], [9, 10]. 

Asan example the assessment of hazard during replacement 
of main engine’s fuel pump is considered. In Fig.3 the func- 
tional variables used for the hazard assessment and the values 

attributed by the system’s user, are presented. 


In the considered case, the following values 
were attributed to the particular functional variables : 


+ V3 (range -50++250) — 150 


& 
E2 
2 
5 
= 


hazard level 


Fig. 4. The resulting fuzzy set and final value 
of hazard level due to functional factors . 


CONCLUSIONS 


On the basis of the obtained investigation results 
the following conclusions can be formulated : 


O The application of fuzzy logic makes it possible to estimate 
operator’s safety already during design process of ship po- 
wer plant, especially in its early stages when information on 
hazardous and noxious factors to operator are not precise 
and associated with high uncertainty 


O Incontrast to the classical solutions the assessment of degree 
of hazard to operator by using fuzzy logic rules models more 
realistically qualitative aspects of human knowledge, and 
the inference process itself does not require performing any 
quantitative analyses 


O The application of the special computer software which 
helps in using fuzzy inference, to a great extent decre- 
ases range of efforts connected with building the base of 
knowledge, as well as its size within the computer-aided 
hazardous zone identification system. 


NOMENCLATURE 


a,b,c — linguistic variables 
iC: — value of functional variable for m-th structural unit 


jm 


c — value of service variable for n-th procedure 

f(x) | — function providing a membership value in [0,1] fuzzy set 

Kam — number of service operations appearing in n-th procedure 
i for m-th structural unit 

W, — weighing factor for n-th procedure 

Ww, — weighing factor for j-th functional variable, 

u&x) —membership function of a given element in A fuzzy set 


A, Al, A2, B, B1, B2, C1,C2 — fuzzy sets 


SA — index of influence of endangering and noxious factors 
? on the creating of potential hazard to operator, which 
occur during carrying out service operations realized 
in compliance with the considered n-th procedure 
concerning m-th structural unit 
Lae cs — index of influence of service factors on the creating of 


potential hazard to operator, which occur during carrying 
out service operations realized in compliance with the 
considered n-th procedure concerning m-th structural unit 
I — index of influence of functional factors on the creating of 
potential hazard to operator, which occur during carrying 
out service operations realized in compliance with the 
considered n-th procedure concerning m-th structural unit 


J — number of functional variables 

M — number of structural units distinguished by the considered 
procedures. 

N — number of considered procedures 

V1+V6 — functional variables 

X — numerical space of consideration 
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An analysis of possible assessment of hazards 
to ship shaft line, resulting from impulse load 


Andrzej Grzadziela 
Polish Naval University 


ABSTRACT 


This paper presents a proposal of identification of a degree of hazard to ship shaft line 
due to impulse load caused by underwater explosion. A theoretical analysis was made 
of influence of changes in co-axiality of shafts on second-kind critical velocities resulting 
from elastic deformations of hull structure in vicinity of shaft bearing foundations. Results 
are presented of pilotage tests of underwater explosion performed on training water area. 
A preliminary mathematical model of underwater explosion is given with taking into ac- 
count mass of explosives and distance from the tested object. The drawn conclusions show 


away of identification of hazard to shaft line by using spectral analysis and time - runs of vibration signals 
recorded at shaft bearings. 


Keywords : ship shaft lines, technical diagnostics, modelling, vibrations, underwater explosion 


INTRODUCTION 


Ship propulsion systems are subjected to specific sea loads 
due to waving and dynamical impacts associated with mission 
of a given ship. Sea waving can be sufficiently exactly mo- 
deled by means of statistical methods. Much more problems 
arise from modelling impacts due to underwater explosion. In 
operation of contemporary technical objects including naval 
ships greater and greater attention is paid to such notions as : 
time of serviceability, repair time, maintenance and diagnosing 
costs [1]. Diagnosing process has become now a standard 
procedure performed during every technical maintenance. Out 
of the above mentioned, the notions of time of serviceability 
and maintenance costs seem to be crucial for the diagnosing 
process of ship power plant. Knowledge of a character of im- 
pulse loading which affects ship shaft line, can make it possible 
to identify potential failures by means of on-line vibration 
measuring systems. This way elimination of costly and time- 
-consuming overhauls on dock leads to lowering operational 
costs and increasing ship fighting merits. 


ANALYSIS OF FORCES ACTING 
ON SHAFT-LINE BEARINGS 


Ship shaft lines are subjected to loads in the form of forces 
and moments which generate bending, torsional and axial vi- 
brations. In most cases strength calculations of driving shafts 
are carried out by using a static method as required by majo- 
rity of ship classification institutions. Moreover they require 
calculations of torsional vibrations which have to comply with 
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permissible values, to be performed. Calculation procedures of 
ship shaft lines generally amount to determination of reduced 
stresses and safety factor related to tensile yield strength of 
material — Fig.1. 


The above mentioned methods do not model real conditions 
of shaft-line operation, which is confirmed by the character 
of ship hull response, i.e. its deformations under dynamic 
loads. Much more reliable would be to relate results of the 
calculations to fatigue strength of material instead of its yield 
strength [5]. 


Yield strength 0.000e 


Fig.1. Simulated static bending stresses 
in propeller shaft due to weight of propeller . 


In static calculation procedures no analysis of dynamic exci- 
tations, except torsional vibrations, is taken into consideration. 
In certain circumstances the adoption of static load criterion 
may be disastrous especially in the case of resonance between 
natural vibration frequencies and those of external forces due 
to dynamic impacts. To analyze the dynamic interaction a sim- 
plified model of shaft line is presented below, Fig.2. 


4 À x 
Fig. 2. A simplified shaft-line model for critical speed calculation [4] . 


Let us note : M, - torque, M, - anti-torque. The system can 
be represented by the following set of equations : 


mh+kh= me(@sing+ $ cos p) 


mv+kh= me(—Gcos@ + sing) (1) 


(J + me”) = me(hsin g- Vcos@) +M; -M, 


The presented form of the equations is non-linear. Con- 
sidering the third of the equations (1) one can observe that 
the variables h, v and ọ are mutually coupled. It means that 
any bending vibration would disturb rotational motion of the 
shaft. The third of the equations (1) can be written also in the 

equivalent form as follows : 


J@=ke(vcosp—hsing)+M,-M, (2) 
To obtain the shaft angular speed Q = è constant 
to use time-variable torque is necessary : 
M=M,-M, = ke(hsin@-—vcos@) (3) 


Theoretical analysis indicates that shaft bending deforma- 
tion continuously accumulates a part of shaft torque. However 
the quantity of torque non-uniformity is rather low since shaft- 
line eccentricity is low; it results from manufacturing tolerance, 
non-homogeneity of material, propeller weight and permissible 
assembling clearances of bearing foundations. The condition 
makes it possible to predict that run-out of propeller shaft may 
also happen at other rotational speeds than the critical 1* kind 
speed calculated during design process. Taking into account the 
torque variability one can express Eq. (1) as follows: 


h+ osh = eQ’ coso 
V+ mph = eQ? sing 
G=7M(t) 


On assumption that the average angular speed of propeller 
shaft equals Q, the angle @ will change with time according 
to the equation : 


(4) 


P= Qt+ , (5) 
and the torque M(t) can be represented as follows : 
M(t) = M,—M)= ag: cos(BOt + ng) (6) 


For ship propulsion system the torque pulsation expressed 
by means of Fourier series is much more complex. It addition- 
ally contains components resulting from number of propeller 
blades, kinematical features of reduction gear as well as distur- 
bances from main engine and neighbouring devices. In general 
case occurrence of only one harmonic does not change reason- 
ing logic. After solving the third of Eqs.(4) and substituting 
Eq.(5) to the obtained expression one obtains the following: 


a 

FOT cos(Qt +) +C, +C, 

If Q = const then C, = 0. Similarly: C, = 0 as C, deter- 
mines displacement of the vibration centre axis but does not 
change character of shaft motion, hence to neglect both the 
integrating constants is possible. Assuming additionally that 
generated vibrations are of low energy as compared with the 
power transmitted by the shaft, i.e. cosn = 1 and sin n ~ n, 
and using trigonometric transformations one can obtain the 
following [4]: 


= (7) 


h+ oh = eQ’cos@ + 


+ T {sin[(B + )Qt + ng] -sin[(B - Qt + ngl} 
T (8) 
, V+ MV = eQ coso + 
z T {cos[(B +DOt + ng] + cos[(B — HOt + ng]} 


From Eqs. (8) one obtains three harmonic excitations : Q, 
(B+ 1)Q and(B — 1)Q. Hence assuming that the torque pulsa- 
tion occurs with the frequency B - Q one obtains three cases of 
resonance for one degree of freedom of propeller shaft: 


Okr = Mp 


1 
OQO kr= pe (9) 


1 
OQkrR= pri” 


where : Q kr > Qpe - critical rotational speeds of 2" kind. The 
above given considerations deal only with one torque pulsa- 
tion harmonic. In practice, at a greater number of excitations 
including instantaneous elastic deformations due to impulse 
loading, one should expect a series of critical speeds of 2"4 
kind whose influence will be reduced by damping due to water, 
stuffing-box and bearing supports [6]. 

For long shaft lines of ships the influence of gravity forces 
on critical speeds should be taken into consideration [7]. Accor- 
ding to Eq. (10) the generated vibrations will be then performed 
respective to static deflection axis of the shaft. 


mh+kh = meQ’cos@ 


mV+kv = meQ’sing (10) 
o=0 
Hence the equations obtain the following form: 
h+ wh =e- Q coso 
V+ oiv =e: Q’sing-mg (11) 


J@=-mgecos@ 


Since in the third equation of the set (9), i.e. that for Qro 
appears the exciting torque of the frequency/angular speed ratio 
B = 1 it means that one has to do with the critical state of 2™ 
kind for B = 1, namely : 


(12) 


Occurrence of such kind vibrations is conditioned by non- 
-zero value of e, which — in the case of ship shaft line — appear 
just after dislocation of a weight along ship, a change of ship 
displacement or even due to sunshine operation on one of ship 


Okro) to 
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sides. A similar situation will happen when e varies due to dyna- 
mic excitations resulting from e.g. sea waving or explosion. In 
this case the critical speed will vary depending on instantaneous 
value of e and damping. 

Theoretical analysis of operational conditions of interme- 
diate and propeller shafts indicates that static and dynamic loads 
appear. In a more detailed analysis of dynamic excitations of 
all kinds the following factors should be additionally taken 
into consideration : 


> disturbances coming from ship propeller (torsional, bending 
and compressive stresses); 

> disturbances from propulsion engine (torsional and com- 
pressive stresses); 

> disturbances from reduction gear (torsional stresses); 

> disturbances from other sources characteristic for a given 
propulsion system or ship mission. 


THEORETICAL BACKGROUND 
OF UNDERWATER EXPLOSION 


Information on potential hazard resulting from underwater 
explosion is crucial not only for ship’s commander during 
warfare but also for ship structure designers. Knowledge of 
loads determined during simulative explosions is helpful in di- 
mensioning ship’s hull scantlings [3]. Another issue is possible 
quantification of explosion energy as well as current potential 
hazard to whole ship and its moving system. 

From the point of view of shock wave impact on shaft line, 
underwater and over-water explosions should be considered in 
two situations : 


> when shock wave (or its component) impacts screw pro- 
peller axially 

> when shock wave (or its component) impacts screw pro- 
peller perpendicularly to its rotation axis. 


The axial shock-wave component affects thrust bearing and 
due to its stepwise character it may completely damage slid- 
ing thrust bearing. Rolling thrust bearings are more resistant 
to stepwise loading hence they are commonly used on naval 
ships [3]. The shock wave component perpendicular to shaft 
rotation axis is much more endangering. 

Shock wave can cause: damage of stern tube, brittle cracks 
in bearing covers and tracks, plastic displacement of shaft 
supporting elements including transmission gear and main 
engine, and even permanent deformation of propeller shaft. 
For that reason contemporary naval ships are tested regarding 
their shock strength - irrespective of their mission — already 
during sea trials, Fig.3. 


Fig. 3. Experimental tests of shock strength 
of frigates against underwater explosion . 
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The problem of influence of sea mine explosion on hull 
structure is complex and belongs to more difficult issues of ship 
dynamics. Underwater explosion is meant as a violent upset of 
balance of a given system due to detonation of explosives in 
water environment. The process is accompanied with emission 
of large quantity of energy within a short time, fast running 
chemical and physical reactions, emission of heat and gas prod- 
ucts. The influence of underwater explosion does not constitute 
a single impulse but a few (2 to 4) large energy pulsations of 
gas bubbles [2, 8, 9]. The pulsation process is repeated several 
times till the instant when the gas bubble surfaces. Hence the 
number of pulsations depends a.o. on immersion depth of the 
explosive charge. The character of changes of pressure values in 
a motionless point of the considered area is shown in Fig.4. 


tp 


p max 


Level of hydrostatic pressure 


Z 


t kl tkol 
<> 


impuls nr 1 


amax 


impulsi 
nr 2 


Fig. 4. Run of changes of shock wave pressure and ship hull acceleration 
measured on hull surface during underwater explosion . 


In the subject-matter literature can be found many formulae 
for determining maximum pressure value, based on results of 
experiments, however data on a character of pulsation and its 
impact on ship structures are lacking. 


EXPERIMENTAL TESTS 


To identify underwater explosion parameters a pilotage test 
was performed with the use of the explosive charge having the 
mass m = 37.5 kg. The schematic diagram of the experiment 
is shown in Fig. 6. During the test were measured vibration 
accelerations of casings of intermediate and thrust bearings 
in the thrust direction and that perpendicular to shaft rotation 
axis. The ship course angle relative to the explosion epicentre 
was 45° and the shaft line rotated with the speed n,„ = 500 
rpm. Ship’s distance from the mine and its immersion depth 
was determined by using a hydro-location station and ROV 
underwater vehicle. The vibration gauges were fixed over the 
reduction gear bearing (p. 1) as well as on the intermediate shaft 
bearing (p. 2) as shown in Fig.5. The measurement directions 
(X and Z axes) are presented in Fig.6. 

The performed test was aimed 
at achieving information dealing with : 
= character of shock wave impact on shaft-line bearings, in 
the form of recorded vibration parameters 
= assessment of time-run of vibration accelerations with 
taking into account dynamic features of the signals in set 
measurement points 


= assessment of possible identification of influence of pul- 2) 
sation of successive gas bubbles during the time-run of 
vibration accelerations 

= identification of features of the signals by means of spectral 
analysis. 


Fig. 5. Schematic diagram of the propulsion system of the ship in question. 
1 - measurement point of 3-dimensional vibration at reduction gear; 


Since the mass of the explosive charge was small, to reliably 
identify the effect of only first and second pulsation was possi- 


GTNT=37, 5kg 
hw=32m 


Rw=420m 


2 - measurement point of 3-dimensional vibration 
at shaft-line intermadiate bearing 


ble during the test. The time lag of the recorded signals was the 


same in all measurement points, as shown in Fig.7 and 8. 


RMS [m/s?] 


a [m/s?] 


12 


Fig. 6. Schematic diagram of the performed experimental test . 
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Explosion, Port side(LB), Thrust bearing, X axis 


Fig. 7. Rms values of vibration accelerations recorded on port side (LB) and starboard (PB) intermediate bearings and thrust bearings 


Explosion, Starboard (PB), Intermediate bearing, Z axis 
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Fig. 8. Time-run of vibration accelerations in two measurement points . 


The most important information was the assessment of the ratios of maximum and minimum values 
of vibration accelerations measured before, during and after underwater explosion. 
The performed analysis concerned all the measurement points, Tab. 1. 
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Tab. 1. Ratios of vibration accelerations in the measurement points K1 through K8 


Indices 1 min before explosion 


Measurement points: 


a. 
max © min 


Measurement points: 
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an. 
max © min 


Measurement points: 


/a 


max min 
Measurement points: 


/a 


max min 


3.60 


Indices 1 min after explosion 
Measurement points: 


1.23 
K4 
ana Amin 1 4 3 4 


a. 
max min 


Measurement points: 


0.83 
K3 
0.92 


0.83 
K8 
1.06 


0.87 
K7 
1.08 


Location 


K1 — reduction gear, 


X axis, PB (port side) 


K2 — reduction gear, 


X axis, LB (starboard) 


K3 - reduction gear, 
Z axis, PB (port side) 


K4 - reduction gear, 
Z axis, LB (starboard) 


Location 


bearing, X axis, 


It can be observed that values of the a, /a_,, ratios from before and after the explosion are close to each other; it means that 


PB (port side) 


K5 — intermediate shaft 


bearing, X axis, 
PB (port side) 


K6 - intermediate shaft 


K7 - intermediate shaft 
bearing, X axis, 
PB (port side) 


K8 - intermediate shaft 
bearing, X axis, 
PB (port side) 


the explosion impulse caused elastic deformations and next the propulsion system comes back to the initial technical state. 


MODELS OF EXCITATION DUE TO UNDERWATER EXPLOSION 


Analysis of dynamic impacts including impulse ones should take into account basic parameters which influence character of 
time-run of a given signal as well as its spectrum. The basic parameters which identify impulse impact resulting from explosion, 


are the following : 


* form of impulse which identifies kind of impulse (Fig. 9a through 9d) 


Excitation time-run At, = 1.36 ms 


Excitation spectrum At, = 1.36 ms 


ae |4 
E | 2E | 
cm i 
|. a ae OE. E, eee 
WAA: 
| flan : 
| | | j \ INS : ARA 
-1 i i i i o. y: VNV V VN ENENININS AEN 
0 0.005 0.01 0.015 0.02 0.025 0.03 0 2000 4000 6000 8000 10000 
t [s] [Hz] 
Fig. 9a. Stepwise impulse. Fig. 9b. Stepwise impulse spectrum. 
l Excitation time-run Atı = 1.36 ms x10 Excitation spectrum At, = 1.36 ms 
: : : 15 r r r r 
AN 
= lt f : \ OE a hgh cate Ms Na tt eh AN este EA EAA ET 4 
= |i | 
& } | 
ai } 
< | } 
0.5 | een | nee nn T es = 
j | 
me 
\ { 
i i i i i VN ee : 
- i 0 ¥ aan, A a cee o —_— 
0 0.005 0.01 0.015 0.02 0.025 0.03 0 2000 4000 6000 8000 10000 
t [s] f [Hz] 


Fig. 9c. Symmetrical impulse. 
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Fig. 9d. Symmetrical impulse spectrum. 


* impulse duration time At, at the ratio A/At 


I I 


maintained constant, which identifies explosive charge power 


(time of propagation of gas bubble) - Fig. 10a through 10d 


Excitation time-run At, = 0.68 ms 


A [mm/s] 


0 0.005 0.01 0.015 0.02 0.025 0.03 


Fig. 10a. Simplified run of shock wave model 
in function of time (At, = 0.68 ms). 


Excitation time-run At, = 0.68 ms 


0 0.005 0.01 0.015 0.02 0.025 0.03 


Fig. 10c. Simplified run of shock wave model 
in function of time (At, = 1.36 ms). 


x10 Excitation spectrum At, = 0.68 ms 


0 i i i T 
0 2000 4000 6000 8000 10000 


Fig. 10b. Spectrum of simplified run of shock wave model 
in function of time (At, = 0.68 ms). 


x10“ Excitation spectrum At; = 0.68 ms 


0 ; : : : 
0 2000 4000 6000 8000 10000 


Fig. 10d. Spectrum of simplified run of shock wave model 
in function of time (At, = 1.36 ms). 


* influence of damping on spectrum form, which identifies distance from explosion 
and - simultaneously - epicentre depth (Fig. 11a through 11d) 


Excitation time-run At; = 0.88 0.4 0.21 ms 


-1 i ; i i i 
0 0.005 0.01 0.015 0.02 0.025 0.03 


Fig. La. Signal which contains transformed impulses 
for amplitude ratios 1: 0.37 : 0.05. 


Excitation time-run At; = 0.88 0.4 0.21 ms 


E E a a 


“0 0.005 0.01 0.015 0.02 0.025 0.03 


Fig. 11c. Signal which contains transformed impulses 
for amplitude ratios 1: 0.18 : 0.02. 


25 Žž 10 Excitation spectrum At;=0.88 0.4 0.21 ms 


0 i l i ; 3 Á i i 
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 


Fig. 11b. Spectrum of the signal which contains transformed impulses 
for amplitude ratios 1: 0.37 : 0.05 . 


x 10 Excitation spectrum At; = 0.88 0.4 0.21 ms 


MMM: o EE ei] 


= 


0 : : : : 4 : $ 
0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 


Fig. 11d. Spectrum of the signal which contains transformed impulses 
for amplitude ratios 1: 0.18 : 0.02. 


POLISH MARITIME RESEARCH, No 3/2007 


OPERATION & ECONOMY 


OPERATION & ECONOMY 


* number of excitation impulses, which informs on distance from explosion, 
combined with explosive charge mass (Fig.12a through 12d) 


Excitation time-run At, = 0.68 ms 


0 0.005 0.01 0.015 0.02 0.025 0.03 


Fig. 12a. Run of shock wave model in function 
of time (two impulses) t = 0.68 ms . 


Excitation time-run At, = 0.68 ms 


0.005 0.01 0.015 


Fig.12c. Run of shock wave model in function 
of time (four impulses) t = 0.68 ms . 


x10” Excitation spectrum At, = 0.68 ms 


1.4 


0 2000 4000 6000 8000 10000 
Fig. 12b. Spectrum of run of shock wave model in function 
of time (two impulses) t = 0.68 ms . 


x10° Excitation spectrum At, = 0.68 ms 


nE EA NOE eL i T d 


0 2000 4000 6000 8000 10000 


Fig. 12d. Spectrum of run of shock wave model in function 
of time (four impulses) t = 0.68 ms . 


* time between successive impulses, which characterizes explosive charge mass (Fig 13a through 13d) 


Excitation time-run At, = 0.68 ms 


-1 i i 
0 0.005 0.01 0.015 0.02 0.025 0.03 


Fig. 13a. Run of shock wave model in function of time with taking 
into account effect of three successive gas bubbles (short time) . 


Excitation time-run At, = 0.68 ms 


0 0.005 0.01 0.015 0.02 0.025 0.03 


Fig. 13c. Run of shock wave model in function of time with taking into ac- 


count effect of three successive gas bubbles (four times longer duration time) . 
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x10” Excitation spectrum At, = 0.68 ms 
1.4 : j 


J2 ee ; EEEE ' EET ; bck gnome eee aapea taaa 


0 2000 4000 6000 8000 10000 
Fig. 13b. Spectrum of shock wave with three successive impacts 
of gas bubbles (short time) . 


-3 


x10 Excitation spectrum At, = 0.68 ms 


0 2000 4000 6000 8000 10000 


Fig. 13d. Spectrum of shock wave with three successive impacts 
of gas bubbles (four times longer duration time) . 


The modelling of run of hydrodynamic pressure changes 
due to underwater explosion, in the direction normal to tangent 
plane at any point of hull, can make it possible to identify the 
quantities which characterize underwater explosion. The formu- 
lae for maximum shock wave pressure and pressure values in 
function of time, presented in the literature [2], do not take into 
consideration all the components whose values are significant 
for assessing degree of hazard to propeller shaft operation. As 
observed in Fig.4 every non-disturbed underwater explosion 
(by shallow water, bed structure, water non-homogeneity etc) 
consists of at least three shock wave impacts resulting from gas 
bubble pulsation. In practice the form of shock wave is limited 
by the following factors : 


* if the distance between explosion epicentre and ship hull 
is small then only one time-extended impulse impact may 
be recorded 


Æ ifthe distance between explosion epicentre and ship hull is 
large then impacts of the second or successive pulsations 
are contained within measurement background and their 
degree of danger to propeller shaft operation is negligible 
(for unlimited water area it depends on explosive charge 
mass and distance). 


The possible recording of measured shock wave pressure 
and accelerations on intermediate and propeller shaft bearings 
enables to identify some explosion parameters hence also 
hazards to power transmission system. Analysing the run of 
underwater shock wave pressure one is able to assume its 
time-dependent function (Fig.14a and 14b) : 


0 20 40 60 80 100 
Fig. 14a. Example of the function form for b = 1.5, c = -0.15 and k = 1 


0 20 40 60 80 100 
Fig. 14b. Example of the function form for b = 1.5, c = -0.15 and k = 10 


A= at" - eX (normalized by a) (13) 


For the assumed mathematical model of the first shock 
wave impulse the run of vibration accelerations recorded on 
ship hull - for the example function given in Fig. 14a - can be 
presented as shown in Fig. 15a and 15b. 


Fig. 15a. Run of the assumed vibration acceleration model 


x 10° 


A [mm/s] 
= 
Nn 


0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 


Fig. 15b. Spectrum of the assumed vibration acceleration model 


FINAL CONCLUSIONS 


It’s common knowledge that failure frequency is the most 
hazardous factor in marine industry, just after aeronautics. 
Dynamic reactions which occur on ships in service at sea are 
rarely able to produce wear sufficient to cause a failure. 


The possible application of an on-line monitoring system 
of vibration parameters of the propulsion system of mine 
hunter makes it possible to perform the typical technical 
diagnostic tests of torque transmission system and to iden- 
tify possible plastic deformations of hull plating as a result 
of underwater explosion. 


The modelling of impulse impact form and next its identi- 
fication makes it possible additionally : 


¢ to identify explosion power by using an analysis of the 
first vibration impulse amplitude and its duration time 


¢ to identify distance from explosion epicentre (hence 
a degree of hazard) by analysing signal’s damping 


¢ to identify a kind of explosion and even characteristic 
features of type of used mine 


è toselect dynamic characteristics of a measuring system 
which has to comply with requirements for typical 
technical diagnostics and for a hazard identification 
system 


¢ to identify elastic or plastic deformation of shaft line by 
using spectral assessment of its characteristic features 
from before and after underwater explosion. 


The presented results of modelling related to the performed 
experimental test do not make it possible - due to strongly 
non-linear character of interactions occurring in sea envi- 
ronment — to assign unambiguously the modelled signal 
features to those of the recorded ones during the real test. 
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O Successive experimental tests will make it possible to verify 
features of the signals assumed for the analysis, to be able 
to build reliable models. 


O The wide range of stochastic dynamic loads acting on 
ships during its life-time makes that in the nearest future 
the application of on-line diagnostic techniques to ship 
propulsion systems, based on analysing vibration signals, 
will constitute an obvious tactical and technical necessity. 


NOMENCLATURE 

M, - driving torque 

M, -anti -torque 

m  — mass 

e — shafts axis displacement 
J — moment of inertia 

h, v — coordinates 

Q, œ — angular velocity 

k — stiffness 
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aa C onference ~ 


SymSO 2006 


27" Symposium 
on Ship Power Plants 


On 16-17 November 2006 27" Symposium on Ship 
Power Plants was held in Szczecin. It was organized by 
Division of Ship Engines and Power Plants, Department of 
Heat Machines and Ship Power Plants, Maritime Techno- 
logy Faculty, Szczecin University of Technology. Scientific 
patronage over the Symposium was taken by : the Marine 
Technology Unit of the Transport Technical Means Sec- 
tion, Transport Committee, Polish Academy of Sciences 
(PAN) and the West- Pomeranian Regional Group of the 
Section on Exploitation Foundations, Machine Building 
Committee, PAN. Arrangement of 27" Symposium on Ship 
Power Plants has been included into the program of cel- 
ebration of 15" Anniversary of the Maritime Technology 
Faculty as well as 60" Anniversary of Szczecin University 
of Technology, commencing in 2006. 

Such symposia have been organized each year by simi- 
lar Divisions, Departments and Institutes which educates in 
designing, construction and operation of ship power plants, 
being organizational units of : Polish Naval University 
of Gdynia, Gdansk University of Technology, Gdynia 
Maritime University, Maritime University of Szczecin 

and Szczecin University of Technology. 

The Symposia are aimed at forming a forum for ex- 
change of scientific technical information and experience 
among scientific workers and persons employed in mari- 
time economy enterprises and institutions in the area of 
the designing, manufacturing and operating of ship power 
plants and their devices, as well as of marine environment 
protection and for tightening relations between science 

and practice. 

The organizers had honour to host — apart from rep- 
resentatives of the above mentioned universities — also 
experts in the field of industry and maritime economy. 
It was representatives of H. Cegielski Works in Poznan, 
MAN Diesel A/S, Wartsila Polska Co Ltd, Alfa-Laval 
Polska Co Ltd, PBP Enamor, Bosch Rexroth Co Ltd, as 

well as Szczecin Nowa Shipyard Co Ltd. 

During two-day proceedings 20 papers were presented 
in 7 sessions and another twenty in poster presentations. 
Themes of the papers were assigned to the following 

groups : Design, Diagnostics and Operation, 
and Environmental Protection and Safety. 


The shares of the particular scientific research centres 
in preparation of the papers were as follows : 


ve Polish Naval University 8 
ye Gdansk University of Technology 9 
ye Szczecin University of Technology 5 
we 
x 


Gdynia Maritime University 5 
Szczecin Maritime University 12 
ve Polish Register of Shipping 1 


The next, 28th Symposium on Ship Power Plants 
will be organized by Gdynia Maritime University 
and held in November 2007. 
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Approximation of the index for assessing 
ship sea-keeping performance on the basis 
of ship design parameters 


Tomasz Cepowski 
Szczecin Maritime University 


ABSTRACT 


This paper presents a new approach which makes it possible to take into account sea- 
keeping qualities of ship in the preliminary stage of its design. The presented concept is 
based on representing ship's behaviour in waves by means of the so called operational 
effectiveness index. Presented values of the index were calculated for a broad range of 
design parameters. On this basis were elaborated analytical functions which approximate 


¢ the index depending on ship design parameters. Also, example approximations of the index 


calculated by using artificial neural networks, are attached. The presented approach may 


find application to ship preliminary design problems as well as in ship service stage to assess sea-keeping 
performance of a ship before its departure to sea. 


Keywords : sea-keeping qualities, sea-keeping performance index, approximation, 
ship design parameters, artificial neural networks, rolling, slamming, green-water shipping 


INTRODUCTION 


In ship design process an optimum solution which satisfies 
assumed economical criteria and technical limitations, is sear- 
ched for. The technical limitations contain a.o. performance of 
ship in rough seas, the called sea-keeping qualities. However 
in the preliminary design stage to take into account the whole 
range of ship sea-keeping qualities is very difficult and rather 
inaccurate when using current calculation methods. It results 
from a few problems the most important of them are the fol- 
lowing : 


in the preliminary design stage values of the parameters 
which significantly influence sea-keeping performance of 
designed ship are not yet known (they result from ship hull 
form and mass distribution which are unknown in that stage 
of designing) 

simple accurate relationships between design parameters 
and sea-keeping qualities of ship are unknown 

inthe preliminary design stage, assessment of sea-keeping 
performance of ship is of a descriptive character, left at 
designer’s discretion and hence imprecise as usual. 


In the papers [6 , 7] was proposed a method of approxi- 
mation of selected sea-keeping qualities based on ship design 
parameters and represented by means of amplitude-phase 
characteristics of ship in regular waves. 

On the basis of such approximations it is possible to deter- 
mine statistical quantities of sea-keeping qualities of ship in ro- 
ugh waves and the operational effectiveness index (acc. [2 , 5]) 


which enables to assess quantitatively sea-keeping performance 
of a given designed ship. The operational effectiveness index 
Ey usually expresses probability of the event that ship response 
in given wave conditions will not exceed an assumed level. 
Hence the index Ey takes values from the interval of 0 to 1. 
The higher the index value the better predicted sea-keeping 
performance. In [8] is presented an example of application 
of the above mentioned index (and of the approximation acc. 
[6 , 7]) to selecting optimum design of ship regarding its sea- 
-keeping performance. 

The approximations presented in [6 , 7] make it possi- 
ble to simply and accurately determine transfer functions 
of selected sea-keeping qualities. However, the assessing 
of the sea-keeping qualities on the basis of the operatio- 
nal effectiveness index is associated with the necessity of 
calculation of their statistical values for ship in irregular 
waves and probability of occurrence of the wave parame- 
ters for which ship’s motions exceed an assumed level. 
The general algorithm of calculation procedure of the index 
Ey , based on the transfer functions of sea-keeping qualities 
is presented in Fig.l. The using of the above mentioned algo- 
rithm for calculation of the operational effectiveness index is 
associated with necessity of carrying out many iterations hence 
the approach can not be applied to multi- criterion optimization 
methods as the target function as its form is then too much 
complex. Moreover, to calculate Ey index value, data on sta- 
tistical distributions of wave parameters in a given sea area or 
shipping route should be known, that excludes the method from 
application to fast assessing a given ship design variant. 
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L 


Approximation of the transfer 
functions of sea-keeping 
qualities Y pX] , X3 ... X,) 


| p 


Calculation of the statistical 
quantities of sea-keeping qualities: 
— variances 
— significant amplitudes 
— occurrence frequencies 


ig 


Probability of occurrence of 
the wave parameters 
Hs ,T 


Calculation of the operational 
effectiveness index Ey , acc. [1, 5, 8] 


Fig. 1. Schematic diagram of calculation of the operational effectiveness 

index of sea-keeping qualities, Er . Notation : X4 , X3 ... X, — ship design 

parameters, Hs — significant wave height, T — characteristic wave period, 
Yin — transfer function of sea-keeping qualities. 


CONCEPT 


For the above mentioned reasons in the investigations in 
question another concept of determining the index Ey was 
taken into consideration. To simplify the method described in 
Fig.1, the approximation of transfer functions of sea-keeping 
qualities was replaced by the approximation of the index Ey 
depending on the ship design parameters X4 , Xz ... X, on 
a given shipping route. The schematic diagram of the concept 

is presented in Fig.2. 


Approximation 
of the operational 
effectiveness index 
E; =f (X, , X ... X,) 


Ship design parameters known 
in the preliminary design stage 
Xi, X -o Xh 


Fig. 2. Schematic diagram of approximation of the sea-keeping qualities 
index Er on a given shipping route, where : X1 , X3 ... Xn — ship design 
parameters, f — searched for approximation function 


In such case the effectiveness index Ey 
on a given shipping route can be approximated 
in accordance with the formula (1) : 


Ey, =f(X, , X2 sss Xp 


where : 


(1) 


Xı , X3 ... Xp — ship design parameters 
f— searched for approximation function. 


On the basis of the assumed concept the approximation 
function f may concern all selected sea-keeping qualities or 
only a given one (e.g. rolling or slamming). Values of the 
index Ey , approximated with the use of the function f , can 
be related to : 


a given shipping route or region of operation 

a given season or time interval 

a selected type of ships or their group 

an assumed load condition (ballast, full load) 
assumed criterion values of sea-keeping qualities. 


HHH HH 


The approximation function f can be determined on the 
basis of the set of model values of the ship design parameters 
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X 1, Xz... Xp and values of the effectiveness index Ey calcula- 
ted by using exact methods. The above mentioned approxima- 
tions can be elaborated by using statistical methods. An example 
of application of the artificial neural networks to determine the 
approximation function f is presented below. 


Application of the artificial 
neural networks to approximate 
the operational effectiveness 
index Er 


The function for approximating the ship operational effective- 
ness index Ep can be determined according to the formula : 


x—>yY (2) 
where : 
X — set of assumed ship design parameters 
Y — setofvalues of the ship operational effectiveness index 
Ey calculated by means of exact methods 
f — searched for analytical function in the form of artificial 


neural network, intended for the approximating of the 
index Er. 


It was preliminarily assumed that the investigations in 
question will concern approximation of the sea-keeping qua- 
lities which have the most detrimental influence on ship’s 
safety, to which — in accordance with [2 , 5] — the following 
ones belong : rolling, slamming, green water shipping onto 
the deck, propeller emerging, vertical accelerations at bow 
and bridg, and pitching. And, successive investigations sho- 
wed that for the considered series of ships in assumed wave 
conditions vertical accelerations and pitching impair ship’s 
safety to a small extent only. Hence only the values of the 
operational effectiveness indices which concern the following 
qualities, were taken into account as the elements of the set 
Y in the equation (2) : 


rolling — (Etyot) 

slamming — (Etslam) 

green water shipping — (ETgreen) 
propeller emerging — (Eqprop) 

all the above specified sea-keeping 
qualities in total — (Ey). 


+++4+ 


And, as the elements of the set X in the equation (2) the 
ship design parameters which significantly influence the above 
mentioned sea-keeping qualities, were taken into account. In 
compliance with the literature sources [1, 3, 5] the following 
ship’s parameters were assumed : 


A L  — length 
4 B — breadth 
4 Cb -block coefficient of immersed 
part of ship hull 
A GMp- initial transverse metacentric height 
A T  —draught. 


The process of searching for the best network 
consisted of the following steps : 


determination of the best network 
structure by using genetic algorithms 
teaching the network 

testing the network 

assessing approximation accuracy 

of the network on the basis of test data. 


VVV WV 


For teaching the neural networks at most 50 % of all the data 
was utilized not to result in over-teaching the networks. 


For the assessing of approximation accuracy errors 2.0 


in teaching and testing, which can be determined — Z: Ae e neen 4 
from Eq. (3), were used : 16 — GRIM 4 
& 4 — CTO-WARES 
2 of 
RMS=,|Ctw En) 3) £12 
n & 1.04 
n 
where : 8 0.8 
RMS — value of error ote 
Erw — model (reference) values used 5 0.4 
in teaching or testing the neural network 0.2 
Er — values calculated with the use 0.0 r = 
of the neural network o0 02 0:4 “ i 1:0 l2 ia 
n — number of records. Pregnengyil ls] 


Fig. 3. Comparison of heave transfer functions : By, = 90°, V = 0 kn ; 
Model values 200 000 dwt tanker of the length L = 310 m, breadth B = 47.17 m, 
draught T = 18.9 m [3] 

As model values the operational effectiveness indices calcu- 
lated for the North Atlantic in winter were assumed. Values of 
the indices were calculated for the series of container carriers 
having the design parameters specified in Tab.1. 

Their sea-keeping qualities were calculated by using exact 
numerical methods with the help of GRIM ! software. In Fig. 
3 — 5 are presented tests of accuracy of the GRIM software as 
compared with the results obtained by recognized scientific 
centres (measurements performed in Wageningen model basin, 
calculations by Delft calculations by the WARES software, 
Ship Design and Research Centre, Gdańsk). The presented 
comparison shows very high accuracy and conformity of the 
calculations performed with the help of GRIM software. 

Values of the operational effectiveness indices were cal- 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
culated in accordance with the algorithm presented in Fig.1 Frequency [1/s] 
complying with [2, 5, 8]. Criterion quantities for sea-keeping 
qualities were assumed according to [2]. Calculations of sec- Fig. 4. Comparison of roll transfer functions : B, = 90°, V = 0 kn ; 
ondary effects of ship oscillation motions were performed for 200 000 dwt tanker of the length L = 310 m, breadth B= 47.17 m, 
the points indicated in Fig. 6. draught T = 18.9 m [3] 


Tab. 1. Model data set : 


* the ship design parameters : L — length, B — breadth, Cb — block coefficient of immersed part of hull, 
GM, — initial transverse metacentric height, T — draught; 
* the operational effectiveness indices : E7—for all sea-keeping qualities, E Tyo — for rolling, E Tyiam — for slamming, 
ETprop. — for propeller emerging, E Tgreen — for green water shipping onto the deck. 


a 
= 


o Wageningen-measurements 
— Delft-calculations 

— GRIM 7 
— CTO-WARES 


N 
Ur 


N 
S 


Roll transfer function [-] 


Variant GMpo [m] | T [m] ETgreen 


ZliSlelælxailalan] alajn] — 


= 
N 


j= 
U 


= 
i 


jæi 
Nn 


p= 
a 


1 The GRIM software was elaborated by the Department of Ocean Engineering and Marine Systems , Faculty of Maritime Technology , 
Szczecin University of Technology 
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= 1.0 
Ss o Wageningen-measurements A(x = 0.02 Lpp ,y=0,z=0.5 T) 
D 0.9 — Delft-calculations “| 
a| Tos — GRIM | g 
=) g — CTO-WARES B (x = Lpp , y = 0.25 B , z = H + 0.45 T) 
S| .2 0.7 
Ee} © 
SEN 
So! 0.5 = X 
2 
2 0.4 0 a 
= ; aD A 
£03 É C(x=Lpp,y=0,z=0) ~ 
so y 
20.2 fi 
a 0. * 
0.1 
0.0 
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Frequency [1/s] 
: Fig. 6. Coordinates of the points at which values of the selected sea-keeping 
Fig. 5. Comparison of pitch transfer functions : By = 180°, qualities were calculated : A — for propeller emerging, B — for green water 
V = 0 kn ; 200 000 dwt tanker of the length L = 310 m, shipping onto bow deck, C — for slamming , where : Lpp — ship length be- 
breadth B = 47.17 m, draught T = 18.9 m [3] tween perpendiculars, B — ship breadth, T — ship draught, H — ship depth. 
Test data 


The testing of approximation accuracy of the function fn was carried out on the basis of the following data : 


= those contained within the range of model data (interpolation) — Tab. 2 
= those from outside the range of model data (extrapolation) — Tab. 3. 


The sea-keeping qualities were calculated by using the methods described in [6, 7], 
whereas values of the indices presented in Tab 2 and 3 — in the same way as for the model data. 


Tab. 2. Test data set for interpolation : 


* the ship design parameters : L — length, B — breadth, Cb — block coefficient of immersed part of hull, 
GM, — initial transverse metacentric height, T — draught; 
* the operational effectiveness indices : Ep — for all sea-keeping qualities, E Tyo — for rolling, ETslam — for slamming, 
Etprop. — for propeller emerging, E Tgreen — for green water shipping onto the deck. 


Variant | L [m] ETron Eqsiam ETgreen 


| 8 ~~ | 232.70 


Tab. 3. Test data set for extrapolation : 


* the ship design parameters : L — length, B — breadth, Cb — block coefficient of immersed part of hull, 
GM, — initial transverse metacentric height, T — draught; 
* the operational effectiveness indices : Er —for all sea-keeping qualities, Eqygy — for rolling, Eq,1qm — for slamming, 
ETprop. — for propeller emerging, E Tgreen — for green water shipping onto the deck. 


Variant | L [m] B [m] Cb [-] Etprop. 
1 114.59 | 19.10 | 0.67 0.87 


2 
3 
4 
5 
6 
7 
8 


Approximation 


The MLP network of 5x1x1 structure (Fig.7), characterized by the statistics described in Tab.4, 
appeared the best for approximating the operational effectiveness index Ey for all sea-keeping qualities. 
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of interpolation and extrapolation. Moreover, the function (4) 


O 


Tab. 4. Statistics of regression issues 
for the neural network approximating the index Er. 
: testin testin 
teaching |. GA 8. 
> O (interpolation) | (extrapolation) 


Io Coefficient of 
correlation R 


Fig. 7. Structure of MLP artificial neural network RMS error 
for approximating the index Er. 


The searched for function approximating the index Ep , elaborated by using the above presented neural network, 
is presented by means of Eq. (4), in an analytical form : 


2.5614 — 0.1378 |+ 4.1667 


0.0088-L - 1.4110 
0.0614-B - 1.4110 

-| | 9.0909-Cb -6.0909 |x [-0.6815 -0.0286 -0.1163 0.3014 2.3751]-0.1378 (4) 
0.4785-GMg-1.1531 


5.5556-T — 1.2203 
l+e 
E, = 


5.5556 


where : 


L — ship length, B — ship breadth, Cb — block coefficient of immersed part of hull, 
GMp — initial transverse metacentric height, T — ship draught. 


As results from Tab. 4 and Fig. 8 and 9, the function de- In Tab. 5 are described the artificial neural networks ap- 
scribed by Eq. (4) shows rather high accuracy both in the range proximating the remaining operational effectiveness indices 


of sea-keeping qualities such as : 


has a very simple structure, that makes it possible to use it + rolling (Eron) 
TO 


in the multi-criterion optimization methods based on genetic 


1 
“Ku 


Operational effectiveness index E, 


Operational effectiveness index E, [-] 


1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 


1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 


+ slamming (Etslam) 
+ green water shipping onto the deck (ETgreen) 
+ propeller emerging (ETprop)- 


algorithms. 


The values of statistical parameters speak for relatively high 
approximation accuracy of the networks in question both in the 
range of interpolation and extrapolation. The neural networks 

are also of a relatively simple structure. 


1 Tab. 5. Types, structure and statistical parameters of the neural networks 
for approximating the operational effectiveness indices : E Tyo — for rolling, 
Etiam —for slamming, E Tprop. ~ for propeller emerging, E Tgreen ~ for green 


water shipping, where : R — correlation, RMS — networks error . 
l 2 3 4—5 6 7 8 Type of Mering 
Variant Parameter Structure Teaching | Inter- laevis 


m calculated by using artificial neural networks m model data network 


polation | polation 


Fig. 8. Interpolation of the operational effectiveness index Er 
for test variants given in Tab. 2. 


SUMMARY 
O In this paper was presented a concept of approximation 
1 2 3 4 5 6 7 8 


process of the operational effectiveness index by making 
use of basic ship design parameters. The index - according 
to [1, 5, 8] - can be applied to assess ship sea-keeping per- 


Fig. 9. Extrapolation of the operational effectiveness index Er formance both in the stage of ship preliminary design and 
for test variants given in Tab. 3. in service. 


Variant 
a calculated by using artificial neural networks m model data 
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OPERATION & ECONOMY 


OPERATION & ECONOMY 


Q 


Example approximations of the index by using the artificial 
neural networks are also attached. The presented approxi- 
mations make it possible to determine the operational ef- 
fectiveness indices for the following sea-keeping qualities : 
rolling, green water shipping onto the deck, slamming, pro- 
peller emerging, and all the above given qualities together, 
on the basis of the following ship design parameters : length, 
breadth, design draught, block coefficient of immersed part 
of hull, initial transverse metacentric height. 


The described approximations can be used to select a ship 
design which is characterized by most favourable sea- 
-keeping qualities on a given shipping route described by 
statistical parameters of waves, at given criterion values for 
sea-keeping performance. 


A way of application of the approximations may be the 
same as that described in [8]. 


The approximations make it possible to calculate the ope- 
rational effectiveness index fairly accurately and in a much 
simpler way than in the iterative methods described in 
[6,7]. 


Owing to the simple form of the functions approximating 
particular indices the approximations in question can be 
used as target functions in multi-criterion optimization 
methods in the preliminary ship design stage. They can be 
also used as a simplified method to assess a given design 
variant from the point of view of sea-keeping performance 
(without any necessity of determining ship response to 
wave action and taking into account occurrence probabi- 
lity distribution of wave parameters on a given sea area). 


The presented investigation concept can be also applied to 
approximation of sea-keeping performance indices, carried 
out for : 


¢ a greater group of ships of different types 
¢ various shipping routes and areas 
¢ various seasons of the year. 


The approximation method in question may also find its 
application to the issues associated with ship operation 
in the phase of voyage planning. In this case ship master 
could be provided with useful information on sea-keeping 
performance of his ship on a given shipping route and in 
a given season, well in advance of the ship’s departure, 
without any necessity of carrying out complex calculations. 


NOMECLATURE 


Etslam T 


Enroll z 
ETgreen 


GM, - 


RMS - 


T 
V 
Xi, 
Ym 


Bw 
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— ship breadth 
— block coefficient of immersed part of ship’s hull 
— operational effectiveness index for all sea-keeping 
qualities 
operational effectiveness index for slamming 
— operational effectiveness index for propeller emerging 
operational effectiveness index for rolling 
— operational effectiveness index for green-water shipping 
initial transverse metacentric height 
— ship depth 
— ship length 
multi-layer perceptron 
Naval Ship Research and Development Center Bethesda 
— correlation coefficient 
error of learning or testing the artificial neural network 
— ship draught 
— ship speed 
X3 ... X,— ship design parameters 
— transfer function of sea-keeping qualities 


— wave heading angle relative to ship 
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of shipbuilding industry, that guarantees independence of 
technical ideas in the very important domain of national 
economy. Implementation of the obtained results of rese- 
arch projects is not limited solely to shipbuilding but also 
covers areas of other maritime industry and municipal 
economy, as well as makes it possible to export the gained 
knowledge and advanced technologies. 

CTO takes part in activities of many international so- 
cieties and organizations, initiated establishment of Polish 
Forum of Maritime Industries and co-operates in running 
Polish Technological Platform of Waterborne Transport. 

R&D projects are realized by CTO in co-operation with 
many domestic economy institutions such as shipbuilding 
and ship-repair yards, ship equipment manufacturing firms, 
design offices as well as technical universities and other 
research centres. 

Since 1998 CTO has taken part in the EU Outline 
Programs and realized projects both on one’s own and in 

cooperation with foreign and domestic partners. 


MARINE ENGINEERING 


Modelling the characteristics of axial 
compressor of variable flow passage geometry, 
working in the gas turbine engine system 


Pawet Wirkowski 
Polish Naval University 


ABSTRACT 


INTRODUCTION 


An important problem of operation of ship gas turbine 
engines! is to know influence of changeable technical state of 
engine on its working parameters. Compressor is a unit of gas 
turbine engine especially sensitive on changes in its technical 
state during operation [5,6]. The more-or-less-contaminated 
atmospheric air flowing into the compressor causes a.o. continu- 
ously changing form of blade passages, increased roughness of 
blade surfaces as well as change of mass of compressor rotor. 
It seriously impacts compressor’s operation stability, changes 
its characteristics as well as performance and efficiency of 
the entire engine. If the compressor design contains a control 
system for setting controllable blades of guide vanes (of initial 
swirl guide vane or/and first stage guide vanes) to make co- 
operation of all engine’s units optimal by continuous correcting 
the compressor characteristics then the disturbances occurring 
in work of the system will result in changes in operation of the 
compressor and engine, of a character similar to those due to 
rotational speed changes or contaminated blade passages of 
compressor. 

Individual parametric features of every engine in service are 
identified by using expensive experimental tests conditioned 
by many constructional and operational limitations. Dynamic 
development of computer technique applied already to the 
design stage of engine’s units has made it possible to use it 
also for numerical simulation of changeable technical state 
processes. Such approach greatly shortens the time necessary 
for diagnostic tests aimed at determining a set of ,,defect- 
symptom” relations, as compared with the time-consuming 
and expensive investigations carried out on real objects. Hav- 
ing at one’s disposal an appropriate computer software one is 
able to elaborate the models simulating operation of engine’s 


' Called also further : ,,engines” 


This paper concerns application of mathematical modelling methods to analyzing gas-dy- 
namic processes in marine gas turbines. Influence of geometry changes in axial compressor 
flow passage on kinematical air flow characteristics, are presented. The elaborated math- 
ematical model will make it possible to realize — in the future — simulative investigations of 
gas-dynamic processes taking place in a compressor fitted with controllable guide vanes. 


Keywords : marine gas turbine engine, compressor, modelling of characteristics. 


units, verified within an allowable range of static and dynamic 
loads. The computer software based on mathematical models 
makes it possible to realize numerical experiments consisting 
in putting-in real variables and hypothetical technical states 
of engine. 


CHARACTERISTICS OF AXIAL 
COMPRESSORS WORKING 
IN ENGINE SYSTEMS 


The universal characteristics of compressor (Fig.1), sho- 
wing the dependence of the compression 7, and efficiency 
nof compressor on the air mass flowing through it, m, and 
the rotational speed n, make it possible to determine the most 
favourable conditions of co-operation of the compressor with 
other units of the engine under assumption that parameters of 
the sucked -in air are of the values complying with the so-called 
ISA standard atmosphere (p,, = 1.013 bar, T,, = 273.15 K, 
 =0%). The characteristics serve for the selection of optimum 
conditions for air flow control and assessment of influence of 
operational factors on parameters of the compressor. 

Fig.2 highlights occurrence of unstable work phenomenon, 
on which a schematic diagram is presented of flow round a blade 
of axial stage rotor under motion with the constant rotational 
speed n, for which a change in the air flow rate m is effected. 
Fig.2a shows the schematic diagram of the flow in the condi- 
tions for which the air flow rate corresponds with optimum 
efficiency of the stage. 

The relative velocity vectors w,, w, are parallel to camber 
line of blade profile, that is conductive to laminar flow through 
blade passages. The lowering of air flow (Fig.2b) as compared 
with calculation conditions at the circumferential velocity u 
maintained constant, makes the axial component of absolute 
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velocity, c, , smaller, that results in the increasing of the inlet 
angle i of air flow onto rotor blades. It is conductive to bound- 
ary layer separation out of concave surfaces of blades and to 


generation of whirl zones. 
S) ni N Nmax 
n2 Nnom 
ny 


Ngy 


Ts 


. n nom 
ni 


m [kg/s] 

Fig. 1. Characteristics of axial compressor of gas turbine engine; 
a-a — stable operation limit, b-b — limit of maximum air flow rates, 
c-c — line of optimum values of compressor efficiency, 

d-d — line of compressor-turbine co-operation . 


A similar phenomenon occurs on convex surface of blade 
(Fig.2c) when air flow rate becomes greater at the circumfer- 
ential velocity maintained constant. Then the air flow rate m 
takes its maximum values. 

At critical values of angles of flow through the produced 
whirl zones which constitute circumferentially displacing 
low-pressure zones, a sudden air flow reversal (pumping) 
towards compressor’s inlet can happen, that results in violent 
flow pulsations transferred onto engine’s structure. The phe- 
nomenon is detrimental and dangerous because of thermal 


and mechanical overload of engine’s structure. The occurring 
vibrations of significant amplitudes may cause fatigue cracks 
in blades [1, 2, 4]. 

In this connection the compressor should be — within the 
range of service rotational speeds — so adjusted as to place the 
line of compressor-network co-operation with certain stable 
work margin [5, 6, 9]. 

During engine’s operation, rotational speed of rotor, air flow 
rate and optimum form of kinematics of air flow through stage 
blade passages, determined by the air flow inlet angle i onto the 
blades, have the greatest impact on compressor’s performance 
and efficiency. The main principle of compressor control during 
changing its rotational speed or flow rate is to maintain values of 
the air flow inlet angles i close to zero. One of the used methods 
for control of axial compressor is to change geometry of its flow 
passage by applying a controllable inlet guide vane or controlla- 
ble guide vanes of a few first stages of compression [4]. 


CONTROLLABLE INLET GUIDE VANE 


The application of controllable blades of inlet guide vane 
and guide vanes of particular compressor’s stages makes it pos- 
sible to simultaneously change inlet angles of flow onto blades 
of rotor rings of the stages by changing the setting angles of 
the blades of the guide vanes during compressor’s rotational 
speed changing. Fig.3 shows the essence of control of blades 
of controllable guide vanes by using single compression stage 
as an example. The situation shown in Fig.3b where velocity 
directions and values are indexed “1”, corresponds with average 
values of rotational speed service range of compressor’s rotor. 
In this case takes place an average angular setting of blades of 
guide vane ring, at which the inlet angle of flow onto blades of 
rotor ring does not generate disturbances in the flow through 
blade passages. In the case of smaller values of rotational 
speed of compressor, i.e. appearance of a smaller value of the 
axial component of absolute velocity, c, , it is necessary to 
decrease the outlet angle of flow from the controllable blades 
of the guide vane ring (Fig.3a) to such an extent as to maintain 
the same value of inlet angle of flow onto rotor blades. The 
analogous situation takes place during operation of compressor 
at greater values of rotational speed of rotor, for which the axial 


Cia>Cia obl 


Fig. 2. Schematic diagram of flow round blades of compressor axial stage ring 
at constant rotational speed of the rotor and changeable angles of air flow inlet. 
a) calculation angle of air flow inlet, b) positive angle of air flow inlet of, 

c) negative angle of air flow inlet . 
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component value of absolute velocity, c,,” , increases. Then it is 
necessary to increase the outlet angle of flow from guide vane 
blades, Fig.3c, to maintain stable work of compressor hence 
a constant value of inlet angle of flow onto rotor blades. 


a) c) 


BE 


Ring 
of controllable 


guide vane 


Ring of blades 
of stage rotor 


Fig. 3. Essence of control of compressor 5 axial stage by changing the 
setting angle of blades of guide vane rings at changeable air flow velocity; 
a) decreased axial velocity, b) calculation axial velocity, 

c) increased axial velocity . 


Application of a control system of geometry of flow passa- 
ges to gas turbine engine of a given design type significantly 
influences unsteady processes. In multi-stage axial compressors 
of compression value exceeding 8 + 10 the design solution 
which ensures stable operation, is to apply a.o. controllable 
blades of inlet guide vane [4]. In Fig.4 are presented compressor 
characteristics of the engine fitted with controllable inlet guide 
vane, on which the lines of rotor acceleration and deceleration 
have been plotted. 

The lines depicted on the characteristics represent the 
compressor’s operation at three angular settings of inlet guide 
vane blades, a,,,,, and the compressor rotational speed, n, , 
maintained constant. 


----- J> 


Combustion 


1 
1 
1 

Il Free 

Propulsion i power 

line ı turbine 

t 


blades of stage 


Limit 
of compressor 
stable operation 


Okw = -10° 
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Fig. 4. Compressor characteristics plotted in the system of relative-value 
parameters — test stand measurements. Change of compressor operation 
range resulting from interaction of controllable inlet guide vane [5]; 
ABB’B”’— points of co-operation 
of the compressor with network in stable states; 

ACC’B”’— points of co-operation 
of the compressor with network during engine acceleration; 
B”D’DA — points of co-operation 
of the compressor with network during engine deceleration . 


From the run of the co-operation lines shown on the 
characteristics it can be concluded that in the engine fitted 
with controllable inlet guide vane, to control output parameters 
of its compressor at constant rotational speed of rotor unit, is 
possible. 

The DR76 and DR77 triple-shaft engines operating with- 
in COGAG propulsion system installed on naval ships of 
,,larantula” class, are equipped with such design solution of 
inlet guide vane. In Fig.5 is presented the schematic diagram 
of such engine with indicated control cuts in its flow part. 


=-=- N 


Fig. 5. Schematic diagram of ship gas turbine engine fitted with separate reversible propulsion turbine [6] . 
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CONTROLLABLE GUIDE VANES 
OF FIRST STAGES OF COMPRESSOR 


In modern solutions of currently manufactured marine en- 
gines their compressors are fitted with controllable inlet guide 
vanes and controllable guide vanes of their first stages in order 
to ensure sufficient stable work margins. Such compressors are 
characterized by high compression values exceeding 20. Of that 
type are the objects being a part of e.g. LM 2500 engines used 
for propelling the frigates of Oliver Hazard Perry class as well 
as TW-3 aircraft engines applied on Mi helicopters. 

A characteristic feature of the sixteen-stage axial com- 
pressor of LM 2500 engine is the possibility of changing of 
setting of the outlet flow angles a, of the blades of inlet guide 
vane and those of its first six stages in function of engine load. 
The solution prevents from unstable work of the compressor 
during fast realization of transient processes from one stable 
state to another. In the case of aircraft version of the engine 
such solution makes it possible to transit from “low gas” state 
to “full load” within only five seconds not going beyond stable 
operation zone. 

In Fig.6 are shown elements of the control system of angular 
setting of blades of particular guide vanes of the compressor 
of LM 2500 engine. 


Actuation Actuation 
ring ring 
segments segments 


Lever arm 


Vane 
actuation 
lever 


Fig. 6. Elements of the angular setting control system of blades of particu- 
lar guide vanes of the LM 2500 engine 5 compressor 


Fig.7 presents the relation between the outlet flow angle, a, 
of the controllable blade rings of guide vanes and the engine 
load represented by the rotational speed of gas generator’s rotor, 
n» reduced to the normal conditions (acc. ISA) at the inlet 
to the engine. The load control lever installed on the engine 
changes its position within the angular range A = 39° + -3°, 
along with changing the gas generator’s rotational speed within 
the range n, = 5000 + 10000 rpm, which results in changes of 
the angular setting of air outlet from blades of inlet guide vane 
and guide vanes of particular compressor’s stages within the 
range shown in Fig.7. 


MATHEMATICAL MODEL 
OF CHARACTERISTICS OF AXIAL 
COMPRESSOR WITH CONTROLLABLE 
BLADES OF INLET GUIDE VANE 


The problems of elaboration of sufficiently accurate ma- 
thematical models of axial compressors are associated with 
range of simplifying assumptions which determine accuracy 
of the numerical modelling of real object [8]. Simulating in- 
vestigations require a.o. to convert the usual graphical form of 
compressor characteristics (Fig. 4) into functional one suitable 
for numerical calculations, i.e. to the following form : 

n, =f (th, Dor akw) 


(1) 
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Fig. 7. Range of changes of the angle of air outlet from blades 
of rings of controllable guide vanes of compressor, a, in function 


of the rotational speed n_„ and the angular setting of load control lever, À, 
of the LM 2500 engine [7,8] . 


Ns =f (ty, Mes Cw) (2) 

The obtaining of an analytical form of the functional 
relations (1) and (2) which model real characteristics of 
compressor, at maintained minimum approximation error, is 
associated with the difficulties resulting from complex form of 
the functions. In the range of low rotational speeds, isodroms 
of the characteristics exhibit moderate slopes which correspond 
with x, ~idem, whereas in the range of high rotational speeds 
they show steep sections corresponding with m,,~ idem. 

Hence it seems that the least squares method and multi- 
dimensional polynomial-based regression [3] can be an effective 
way for determining an analytical description of axial compressor 
operation, which guarantees that deviations of the model from 
reality would be contained within limits of measurement error. 

The overall model of compressor is searched for by means 
of the set of regression equations which approximate its uni- 
versal characteristics : 


e e 2 2 
m,=ag+a,;m,,+a,(m,) +a,n,,+a,(n,,) + 


: 2 
+as5mMm,, Ny F+aAgAaK wt a7(Oxw) F (3) 
+agM,,lk wt aN, Ax w 
e . 2 2 
T,= bo +b, m,,+ b,(m,,) +b3n,,+b,(n,,) + 
: 2 
+b;m,„ n, + bpd wt b- (Akw + (4) 


+ bm, Akw+t bon,,.AKw 


Values of the regression coefficients a, , b, are determined 
on the basis of the Gauss-Markov theorem when searching for 
minimum values of the functionals: 


J (49,4),47,43,44,85,4¢,47,4g,a9) = 
S 


z = g (5) 
=), (T4.- Ti.) 
k=l 
J, o,b;,b5,b3,b4,b5,b¢,b7, bg, by) = 
(6) 


=») (ns ~ Nck j 
k=1 


They constitute sums of squares of deviations of the values 
obtained by using the model 7, , N from real values T, » 1s. 
Adequacy of matching the mathematical description of 
characteristics of a given compressor and its real run can be 
assessed on the basis of the following factors : 


a) value of the multi-dimensional correlation coefficient R 
expressed as follows : 


N A n 
I = Y Ya — Yn) 
R=——= - - (7) 
Som- Ye- 
n=l n=l 


F : ~? 
b) value of the remainder variance o“ expressed as follows : 


N 
x2 1 
— 8 
ae ae 2% Y,) (8) 


T Value of R should be as large as possible, and value of 
©“ - as small as possible. It then speaks for an insignificant 
deviation of the model from reality. 

The mathematical model of compressor can be improved 
by adding new terms or replacing the existing ones in the 
equations (3) and (4). Solution of the mathematical model of 
compressor consists in determining values of the coefficients 
in the elaborated regression equations. 

The LP compressor of DR76 engine was used to confirm 
usefulness of the least squares method and multi-dimensional 
polynomial-based regression method for building the mathe- 
matical model of characteristics of the compressor with change- 
able flow passage geometry. In Fig.8 a part of the compressor 
model characteristics is presented for three angular settings 
of controllable inlet guide vane : agw = -10°, agw = 0°, ay, = 


15 16 17 18 19 20 21 22 23 
Reduced mass air flow rate m, [kg/s] 

Fig. 8. Compressor model characteristics for changeable angles of setting 
of controllable blades of initial whirl guide vane; stable operation limit line 
Jor the angle of setting of controllable guide vane, app = -10°, 
stable operation limit line for the angle of setting of controllable 

guide vane , G4, = 0°, stable operation limit line for the angle 

of setting of controllable guide vane , a, = +10°. 


= +10°. Generally, for each angular position of guide vane 
blades the compressor characteristics take the form described by 
three values : of the rotational speed of compressor rotor, n, of 
the mass air flow rate Mp and the compressor compression 7. 
In Fig.8 can be observed changes of the compressor operation 
range for two values of the reduced rotational speed of 
compressor, n, = 8000 and 9000 rpm, in function of the angle 
of setting of controllable blades of initial whirl guide vane. 

Values of the regression coefficients of the equations which 
approximate characteristics of the compressor in question are 
given in Tab. 1, and in Tab. 2 — values of statistical parameters 
showing degree of adequacy of values obtained from model 
investigations against experimental ones. 


Tab. 1. Values of regression coefficients for particular equations. 


Regression Regression 
Ne coefficient Valte coefficient Velie 
1 a, 2.099955 b, -0.583216 
2 a, -0.312917 b, -0.272626 
3 a, -0.022847 b, -0.017450 
4 a, 0.000047 b, 0.000976 
5 a, 0.000000 b, 0.000000 
6 a, 0.000115 b, 0.000113 
7 a, 0.059954 b; 0.068760 
8 a, -0.003556 b, 0.000929 
9 a, 0.016743 b, 0.009397 
10 a, -0.000038 b, -0.000032 
Tab. 2. Degree of adequacy between the model and real object 
No. Characteristics R 5? 
1 n =f (ma: Ms OK w) 0.99939 | 0.00628 
2 | nsf (ma Dies aw) 0.99782 | 0.00298 
CONCLUSIONS 


O The presented mathematical model of axial compressor 
working in the gas turbine engine system makes it possible 
to numerically investigate the gas-dynamic processes taking 
place in its flow passages. 


O The changeable setting angles of blades of controllable 
guide vanes of particular compression stages, taken into 
account in the model, significantly widens possible iden- 
tification of unserviceability states in the system of flow 
passage geometry control. 


O Simulative investigations of gas-dynamic processes taking 
place in the compressor in the conditions of introduced 
changes in set-up values of the input quantities of the pre- 
sented model, make it possible to determine the diagnostic 
relations ,,defect-symptom” used for assessing technical 
states of gas-turbine engines. 


NOMENCLATURE 

a, - angle of air flow off blades of guide vane ring 

eg - setting angle of controllable blades of guide vanes 

Bo B, - angle of air inlet to and outlet off rotor blades, 
respectively 

c - axial component of absolute flow velocity 
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Cia - axial component of absolute flow velocity at inlet to 
rotor blade ring, 

Cii - computational value of axial component of absolute 
flow velocity at inlet to rotor blade ring, 

n, - compressor efficiency 

i - angle of air flow onto rotor blades 

0) - relative air humidity 

À - angular setting of engine load control handle 

m - mass air flow rate 

Mpal - mass fuel flow rate 

n - rotational speed of compressor rotor 

Ny, - rotational speed of rotor of compressor during engine 


idle running 
n, n,n, - rotational speeds of compressor rotor at partial load 


states 

Dom? max 7 FOtational speeds of compressor rotor at rated and 
maximum load of engine, respectively 

n, - reduced rotational speeds of compressor rotor 

Po - ambient pressure 

TT. - compression of compressor 

R - multi-dimensional correlation coefficient 

Ta - ambient temperature 

~a? . . 

om - remainder variance 

u - circumferential velocity 

WW, - relative air flow velocity at inlet to and outlet from rotor 


blade ring, respectively 


Aw - flow swirl in rotor 
Abbreviations and indices: 

a, b, - coefficients of regression equation (i = 0...9) 

COGAG - Combination Gas-turbine and Gas-turbine propulsion 
system 

k - number of regression function coefficient 

KS - combustion chamber 

MN - space between propulsion turbine and LP turbine of 
engine 

n - number of measurement (n = 1 ... N) 

N - set of measurement points 

PM - space between compressors 

PR - reduction transmission gear 

PSK - space between HP compressor and combustion chamber 

S - compressor 

SN - screw propeller 

SNC - low pressure (LP) compressor 

SWC -high pressure (HP) compressor 

TN - propulsion turbine 

TNC - low pressure (LP) turbine 

TWC -high pressure (HP) turbine 

Y, - real value at the instant of n-th measurement 

Yn - value calculated from model at the instant of n-th 


measurement 


Y. - real arithmetic mean value 
Ya - arithmetic mean value calculated from model 
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Miscellanea 


Days of Engineering 


On 23-25 November 2006, 37" Days of Engineering 
was held on the occasion of 60" Anniversary of Technical 
University of Szczecin, the oldest technical university of 
West Pomerania, as well as of many scientific technical so- 
cieties of the town, namely : that of electricians, geodesists, 
mechanical engineers and technicians, water engineers and 
technicians as well as the Federation of Scientific Technical 
Societies NOT (Naczelna Organizacja Techniczna). 


The jubilee was celebrated under the slogan : 


Youth and Engineering — a chance 
to developing the Town and Region 


which has had a very distinct meaning. 
It was arranged by the following organizations : 


the Federation of Scientific Technical Societies NOT, 
Szczecin 

Szczecin Division of the Polish Society of Electri- 
cians 

Technical University of Szczecin 

The Society of Graduates from Technical University 
of Szczecin. 
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